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We study the three nucleon force in the triton channel using dynamical clover fermion lattice 
QCD. The Nambu-Bethe-Salpeter wave function is utilized to obtain the potentials among three 
nucleons. Since the straightforward calculation is prohibitively expensive, two different frame- 
works are developed to meet the challenge. In the first method, we study the effective two nucleon 
potentials in the three nucleon system, where the differences between the effective two nucleon 
potentials and the genuine two nucleon potentials correspond to the three nucleon system effect, 
part of which is originated from the three nucleon force. The calculation is performed using 
Nf = 2 clover fermion at m K = 1.13 GeV generated by CP-PACS Collaboration, and Nf = 2+1 
clover fermion at m K = 0.70,0.57 GeV generated by PACS-CS Collaboration. In the second 
method, we study the three nucleon system with 3D-configuration of nucleons fixed. This en- 
ables us to extract the three nucleon force directly, if both of parity-even and parity-odd two 
nucleon potentials are provided. Since parity-odd two nucleon potentials are not available in 
lattice QCD at this moment, we propose a new general procedure to identify the three nucleon 
force using only parity-even two nucleon potentials. The calculation are performed with Nf = 2 
clover fermion at m K = 1.13 GeV generated by CP-PACS Collaboration, employing the linear 
setup for the 3D-configuration. Preliminary results for the scalar/isoscalar three nucleon force are 
presented. 
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1. Introduction 

Since the celebrated work by Yukawa 75 years ago, the nuclear force has been one of the most 
essential quantities in the development of nuclear physics, and so-called "realistic nuclear poten- 
tials" between two nucleons (2N) are available to date. However, recent precise calculations of 
few-nucleon systems clearly point that the 2N force alone is insufficient to understand the nuclei, 
which calls for three (and/or more) nucleon forces. Actually, the three nucleon force (TNF) is pro- 
jected to play an important and nontrivial role in various phenomena in nuclear and astro physics. 
For the binding energies of light nuclei, attractive TNF is required to reproduce the experimental 
data. On the other hand, repulsive TNF is necessary to reproduce the empirical saturation density 
of symmetric nuclear matter. For the EoS of asymmetric nuclear matter, repulsive TNF is required 
to explain the observed maximum neutron star mass. Recently, it is argued that TNF is responsible 
for the anomaly in the drip line and the nontrivial magic number of neutron-rich nuclei 

Although the experimental/theoretical scrutiny of three nucleon scattering are shedding light 
on the natures of TNF [Q], our knowledge on TNF is still quite limited. Pioneered by Fujita- 
Miyazawa [Q], TNF have been mainly studied from the two-pion exchange picture with the A- 
excitation. In addition, repulsive TNF is often introduced phenomenologically [Qj. Recently, the 
TNF based on chiral EFT is developing but the unknown low-energy constants can be obtained 
only by the fitting to the experimental data. Since TNF is originated by the fact that the nucleon is 
not a fundamental particle, it is essential to study TNF from the fundamental DoF, i.e., quarks and 
gluons. In this proceeding, we report such first-principle calculations of TNF using lattice QCD. 

In the lattice QCD calculation of nuclear forces, it is recently proposed [||] to use the Nambu- 
Bethe-Salpeter (NBS) wave function so that the potential is faithful to the phase shift by construc- 
tion. The obtained 2N potentials are found to have desirable features, such as attractive well at long 
and medium distances, and the central repulsive core at short distance [|], pi]. The method has been 
successfully extended to the hyperon-nucleon (YN) and hyperon-hyperon (YY) interactions [pi |9|]. 

We extend this methodology to the three nucleon (3N) system, namely, in the triton channel. 
Due to the significantly enlarged DoF, the straightforward calculation is impossible. We explore 
two different methods to overcome this problem, one is the study of effective 2N potentials in 3N 
system, and the other is the study of 3N system with fixed 3D-configuration by the linear setup. 

2. Formulation for the effective 2N potential 

Since the detailed formulation to study a 2N system is given in Ref. [JJJ, we discuss the exten- 
sion to a 3N system here. We first consider the equal-time NBS wave function \\r(f,p), which can 
be obtained by the calculation of six point correlation function, 



G aS $ % a>f}'y{?3,t-tQ) = (N a (x u t)Np(x 2 ,t)N r (x 3 ,t) (N' al (t Q )NUt )N'Jt Q )) ), (2.1) 



where r = x\ — X2, p = X3 — (xi +X2)/2 are the Jacobi coordinates, and N denotes either of p or n. 

At the leading order of the velocity expansion of the potentials, the NBS wave function can be 
converted to the potentials through the following Schrodinger equation, 
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where Vxtj^j^ij) with r i; - = Xj —xj denotes the potential between (j, j)-pair, Vtnf(?,P) the TNF, 
\L r = Hp = 2ra;v/3 the reduced masses. If we can calculate y/~(r,p) for all f, p, and if all 

V2N,ij{rij) are available by (separate) lattice calculations for the genuine 2N system, we can extract 
Vtnf(7,P)- Unfortunately, this is not the case: Since both r and p have L? DoF, the calculation 
cost is more expensive by a factor of L? compared to the 2N system. Furthermore, the number of 
diagrams to be calculated in the Wick contraction tends to diverge with a factor of N u !xA^! (N u j 
are numbers of u,d quarks in the system). We also note that not all 2N potentials are available in 
lattice QCD at this moment: Only parity-even 2N potentials have been obtained so far. 

In order to avoid these problems, we consider the effective 2N potential in the 3N system. 
More specifically, we take the summation over the location of the spectator nucleon N(xj), 



and define the effective potential between N(xi) and N(x2) via the effective Schrodinger equation, 



In this calculation, the DoF of p is integrated out beforehand, and thus the calculation cost is 
reduced by a factor of ~ 1/L 3 , compared to the straightforward calculation. Yet, the calculation re- 
mains quite expensive due to the large numbers of Wick contractions, and we use several techniques 
to reduce the calculation cost, e.g., we take advantage of symmetries (such as isospin symmetry), 
and we employ the non-relativistic limit for the source nucleon operator. Similar techniques are 
(independently) developed in the calculations of binding energies of triton/helium nuclei JTc|]. 

As the 3N system, we study the triton channel, / = 1/2, J p = l/2 + . Because the spectator 
nucleon is projected to the S-wave, the possible quantum numbers between the (effective) 2N are 
only 2S+l Lj = l S , 3 S U 3 D U and we can obtain the effective 2N potentials V e ff(r) in parity-even 
channel, i.e., the central V^~ l ^~°, V^~°^j~ l and the tensor Vj~^~ l potentials. We calculate all 
counterparts of these parity-even potentials in the genuine 2N system, V2/v(r), and compare them 
with V e ff{r) to extract the effect of the 3N system. In one sense, V e ff(r) — VzN{r) can be considered 
to be the "finite density effect" in the 3N system. Some of this effect are attributed to the genuine 
2N potential with the nontrivial 3N correlation, and the others are originated by the genuine TNF. 
In this way, we can (indirectly) access the effect of TNF. 

3. Lattice setup and the results for the effective 2N potential 

We employ Nf = 2 dynamical configurations with mean field improved clover fermion and 
RG-improved gauge action generated by CP-PACS Collaboration flTT|]. We use 598 configurations 
at j8 = 1.95 with the lattice size of L? x T = 16 3 x 32, which corresponds to (2.5fm) 3 box in 
physical spacial size with the lattice spacing of a -1 = 1.269GeV. We calculate at the hopping 
parameter of u, d quarks K u d = 0.13750, which corresponds to m n = 1.13 GeV, = 2.15 GeV. 
We use the wall quark source with Coulomb gauge fixing. In order to enhance the statistics, we 
perform the calculation for sources on 16 time slices for each configuration. 



<H?)=Lv^p*)=I>(? 5 p*) 



(2.3) 



J-V 2 r +V eff (r) 0(r)=£0(r). 
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Figure 1: Effective 2N potentials, where 
red, blue, brown points correspond to V C ~A ~ , 
^Ceff 1, ^Teff 1 potential, respectively. 
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Figure 3: Same as Fig. g, but for V c 
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Figure 2: The difference between the effective 
2N and the genuine 2N for V^ 1,5- ' potential. 
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Figure 4: Same as Fig. 0, but for V- 
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In Fig. m we show preliminary results for V e //(r) in the triton channel at t — to = 8. Here, 
the constant shift by energy is not included for the central potentials. What is noteworthy is that 
V e ff(r) are obtained with good precision. This is quite nontrivial, since the S/N usually gets worse 
for more quarks in the system. In Figs. ^^JQ, we plot V e ff(r) — V2n{t) f° r eacn potential. One can 
see that the discrepancy is consistent with zero within error-bar. In particular, the tensor potential 
is best constrained within several MeV statistical error, and there is no indication of the TNF effect. 

One of the possible explanations why the TNF effect is not observed is that the TNF is sup- 
pressed at heavy quark mass. Actually, pion exchange is expected to be strongly suppressed with 
m K = 1.13 GeV in this calculation. Therefore, we investigate the quark mass dependence using the 
configurations with smaller quark masses, generated by PACS-CS Collaboration [|T2|]: Nf = 2+1 
nonperturbatively &{a) improved clover fermion, V = 32 3 x 64, /3 = 1.90, aT x = 2.18GeV. We 
use 399 configurations at K ud = 0.13700, K s = 0.13640 (m„ = 0.70 GeV, m N = 1.58 GeV) and 400 
configurations at K ud = 0.13727, K s = 0.13640 (m % = 0.57 GeV, m N = 1.41 GeV). For both quark 
mass setups, we perform 4 source time slice measurements for each configuration. 

We obtain basically similar results to previous results: For all potentials at either quark mass, 
the differences between V e ff(r) and V2jv(r) are found to be consistent with zero within error-bar. In 
particular, the statistical error of V e //(r) — Vwif) in tensor potential remains less than 5 MeV even 
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for the lightest quark mass. This may indicate the necessity to decrease the quark mass further. 

Another possible explanation is that the TNF effect is obscured by the summation over the 
location of the spectator nucleon. In fact, while the TNF effect is expected to be enhanced when 
all three nucleons are close to each other, such 3D-spacial configurations make small contributions 
in the spectator summation. In order to assess this possibility, we study the 3N system with (at 
closely) fixed 3D-configuration, which is the subject of the next section. 

4. Formulation for the fixed 3D-configuration with the linear setup 



We consider Eq. (2.2) with fixed 3D-configuration of r, p. The advantage of this method is 
that the subtraction procedure of the genuine 2N potentials Vin becomes clearer in principle, and 
thus the TNF can be extracted directly. Furthermore, by choosing small |r|, |p| configurations, the 
TNF effect is expected to be enhanced. The disadvantage of this calculation is the calculation cost: 
It is more expensive by a factor of ^(10)-^(10 2 ), compared to the effective 2N potential study. 
Therefore, we restrict the calculation to limited 3D-conngurations. In addition, efforts are taken to 
speed up the calculation code. As the fixed 3D-configuration, we take the linear setup with p = 0. 
Just for the sake of convenience, we redefine as r — > 2r hereafter. In other words, three nucleons 
are aligned linearly with equal spacings ofr= \r\ in this linear setup. 

The advantage of the linear setup is it's simplicity. Because of p = 0, the third nucleon is 
attached to (l,2)-nucleon pair with only S-wave. Considering the total 3N quantum numbers of 
1 = 1 /2,J F = 1 /2 + , the wave function can be completely spanned by only three bases, which can 
be labeled by the quantum numbers of (l,2)-pair as ^0, Si, D\. Therefore, the Schrodinger 
equation can be simplified to the 3x3 coupled channel equations with the bases of s , i//3 S , 
\j/3 Dl . The reduction of the dimension of bases is expected to improve the S/N as well. 

Unfortunately, even by the calculation of the fixed 3D-configuration (including the linear 
setup), the subtraction of V2N remains nontrivial. As was noted in Sec. ||, the parity-odd poten- 
tials are not available in lattice QCD at this moment, and we cannot subtract them unambiguously. 
Note that although the total parity of the 3N system can be projected, a 2N-pair could be either of 
positive or negative parity. The familiar procedure of partial wave expansion cannot be performed 
here, since we can calculate only limited 3D-conngurations due to the huge calculation cost. 

In the effort to overcome this issue, we find that the following channel in the triton is useful, 



p-n n + p n n- - 1111 p- + n n p- + n p n - n p-n- 



(4.1) 



This wave function itself has been well known, but the point here is that it is anti-symmetric in 
spin/isospin spaces for any 2N-pair. Combined with the Pauli-principle, it is automatically guaran- 
teed that any 2N-pair couples with even parity only. Therefore, we can extract the TNF unambigu- 
ously in this channel, without the information of parity-odd 2N potentials. Note that no assumption 
on the choice of 3D-configuration is imposed in this argument, and we can take advantage of this 
feature for the future TNF calculations with 3D-configurations other than the linear setup. 

Coming back to the linear setup, we examine the explicit form of the potential matrix of V%n- 
At the leading order of the velocity expansion, V2N can be written in terms of center V^ s and tensor 
V^ s potentials with isospin / and spin 5, V°°, V^ , V" 1 , V^ 1 , V$ x , V} 1 , where the label "2N" is 
omitted for simplicity. Explicit calculation gives us 
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(4.2) 



Here, we span the spaces with the rotated bases given by (1//5, 1/%, Y 3 Di ) T ' where 1//5 in Eq. ( [Q| ) is 
shown to be 1//5 = ^ ( — So + i//3 Sl ) , and 1/% = -j= ( + So + i/fi Sl ) . Note that neither of parity-odd 
2N potentials, V^°, Vq , Vj- , appear in the first row in Eq. (4^2), as was discussed previously. 



5. The lattice QCD results for the linear setup 

We employ the CP-PACS Nf = 2 clover fermion configurations [O], which are used for the ef- 
fective 2N potential study in Sec. ^. 598 configurations are used with 16 source time slice measure- 
ments for each configuration. We perform the calculation with the linear setup at 7 physical points 



of the distance r. As is explicitly shown in Eq. (|4.2|), we have only one channel which is free from 
parity-odd 2N potentials. Correspondingly, we can determine one type of TNF. In this proceeding, 
we consider the scalar/isoscalar type TNF. In fact, in the Urbana IX model [Q], which is an often- 
used phenomenological TNF, TNF consists of two parts, one is the two-pion exchange TNF (which 
has spin/isospin dependencies), and the other is phenomenologically introduced scalar/isoscalar re- 
pulsive TNF. Since the two-pion exchange TNF is expected to be suppressed in the current lattice 
setup of m % = 1.13 GeV, it is reasonable to assume the scalar/isoscalar nature for TNF here. 

In Fig. |5| we plot each wave function of 1//5, Ym^ Y 3 D[ i n tne triton channel at t — to = 8. 
We observe that \ff s dominates the wave function. This is because \f/ s contains the component for 
which all three nucleons are in S-wave. We emphasize that the observation of good S/N for the 
wave function is much more nontrivial than the effective 2N study, because the practical statistical 
sampling number becomes much smaller by fixing the 3D-connguration. 

By subtracting the V2N in Eq. ( |4~2| ) from the total potentials in the 3N system, we determine the 
TNF. In Fig. |6[ we plot the preliminary results for the scalar/isoscalar TNF. Here, the r-independent 
shift by energies is not included, and thus about <??(10) MeV systematic error is understood. There 
are various physical implications in Fig. |6[ At the long distance region of r, the TNF is small as is 
expected. At the short distance region, we observe the indication of repulsive TNF. Recalling that 
the repulsive short-range TNF is phenomenologically required to explain the saturation density 
of nuclear matter, etc., this is very encouraging result. Of course, we note that further study is 
necessary to confirm this result, e.g., the study of the ground state saturation, the evaluation of the 
constant shift by energies, the examination of the discretization error. 

6. Summary 

We have studied the three nucleon force (TNF) in the triton channel in lattice QCD, developing 
two different methods. In the first method, effective 2N potentials have been studied using Nf = 2 
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Figure 5: The wave function with linear setup 
in the triton channel. Red, blue, brown points 
correspond to y/j, \j/M, V^d > respectively. 



Figure 6: The scalar/isoscalar TNF in the tri- 
ton channel, plotted against the distance r in the 
linear setup. 



clover fermion at m K = 1.13 GeV, and Nf = 2+1 clover fermion at m % = 0.70,0.57 GeV. The 
effective 2N potentials have been found to be consistent with genuine 2N potentials within error. 
In the second method, we have fixed the 3D-configuration of three nucleons. In particular, we have 
established the general procedure which can identify the TNF without the information of parity- 
odd 2N potentials. The calculation have been performed with Nf = 2 clover fermion at m n = 1.13 
GeV with the linear setup for the 3D-configuration, and the indication of repulsive TNF at short 
distances have been obtained. Further work to confirm the results is currently underway. 
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